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ERROR  SiTIB':'  OF  A  WORLDWIDE  SATELLITE  TRIAIJGULATIOR  NET 


ABSTRACr 

A  proposed  three-dimensional  geodetic  reference  system,  consist¬ 
ing  of  36  stations  so  situated  on  the  surface  of  the  earth  as  to  form 
a  closed  system  of  approximate  eq^’uilateral  triangles  surrounding  the 
globe,  is  studied  in  an  attempt  to  find  out  to  vhat  degree  of  accuracy 
the  station  positions  can  be  determined  vhen  triangulating  the  corres¬ 
ponding  station  coordinates  by  photogi-ainmetric  techniques.  The 
polyhedron  of  the  spatially  oriented  triangles  is  formed  by  photo¬ 
graphing  a  high  altitude  satellite  against  the  star  backgroimd.  The 
triangulation  method  depends  on  applying  analytical  photogrammetric 
techniques  to  satellite  observations  where  the  satellite  serves  as  a 
high  altitude  auxiliary  triangulation  target  point.  Numerical 
analyses  on  the  whole  net  and  portions  of  it  are  presented  from 
whicii  the  problem  of  error  propagation  and  the  influence  of  the 
number  of  baselines  can  be  studied . 
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I .  INTRODUCTION 


The  use  of  artificial  satellites  in  scientific  research  provides 

geodesy  for  the  first  time  in  its  history  with  a  means  to  create  a 

three-dimensional  worldwide  reference  system  without  referring  to  the 

direction  of  gravity.  Many  of  the  basic  difficulties  encountered  in 

a  classic  geodetic  triangulation  will  be  eliminated  when  photogram- 

metric  methods  are  applied  to  satellite  observations.  A  satellite 

triangulation  can  be  used  for  establishing  a  worldwide  three- 

dimensional  geodetic  reference  system  which,  without  referring  to 

1* 

gravity,  is  entirely  geometrically  defined.  The  system  will  consist 

of  stations  on  the  physical  surface  of  the  earth,  located  in  such  a 

way  as  to  form  a  closed  system  of  as  near  as  possible  equilateral 

triangles  encompassing  the  globe.  These  stations  will  serve  as  a 

reference  frame  to  which  can  be  tied  both  geometrically-oriented 

mapping  programs  and  the  evaluation  of  satellite  orbits  for  determin- 

2 

ing  gravitational  and  related  geo-physical  parameters. 

Coast  and  Geodetic  Survey  gave  considerable  thought  to  the 

requirements  in  planning  a  worldwide  triangulation  system.  This 

organization  has  the  technical  responsibility  for  the  execution  of  a 

worldwide  satellite  triangulation  program.  Their  conclusions  concern- 

■3_k 

ing  a  usefxil  satellite  orbit  are  summarized  below. 

First,  the  line  of  sight  to  the  satellite  roust  be  at  least  30° 
above  the  horizon  in  order  to  avoid  disturbing  refraction  anomalies. 

Second,  the  maximum  angle  between  any  two  planes  generating  the 
spatial  direction  between  two  stations  must  be  at  least  60°  in  order 
to  provide  sufficient  geometrical  strength  for  the  determination  of 
three-dimensional  station  positions . 

Third,  the  distance  to  the  satellite  should  be  about  the  same  as 
that  between  any  tv/o  stations,  so  that  the  scale  does  not  contribute 
an  unproportional  amount  of  error. 


*  Superscript  numbers  denote  references  which  may  he  found  on  page  72. 
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A  nominal  polar  orbit  is  necessary  because  the  satellite  must 
pass  over  every  geographical  position.  Furthermore^  a  nominal  circular 
orbit  is  required  in  order  that  cheuiges  in  altitude  are  minimized.  A 
near  optimijin  scheme  was  found  which  consists  of  thirty- six  stations. 
Table  I  is  a  list  of  the  names  and  geographical  locations  of  the 
thirty-six  stations,  which  make  up  the  proposed  worldwide  satellite 
triangulation  net.  Figures  1,  2  and  3  show  details  of  the  thirty-six 
station  net  arrangement  on  the  globe. 

An  error  analysis  of  the  proposed  thirty- six  station  net  is 
conducted  to  determine  with  what  accuracy  the  station  positions  can 
be  obtained  by  satellite  triangulation.  The  reference  system  with 
respect  to  which  individual  satellite  directions  are  determined  is 
the  right  ascension  and  declination  system.  The  triangulation  is 
accomplished  by  photographing  a  high  altitude  satellite  against  the 
star  background  simultaneously  from  at  least  two  ground  stations . 
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TABLE  I 


PROPOSED  V70RLDWIDE  SATELLITE  TRIANGULATION  NET 


No. 


Station  Name 


Latitude  Longitude 


1 

Greenleuid,  Thule  A.F.B. 

76.5° 

N 

68.7° 

w 

2 

U.S.A.,  Aberdeen,  Hd. 

39.5 

N 

76.1 

\] 

3 

U.S.A.,  Larson  A.F.B.,  Wash. 

47.2 

N 

119.3 

\1 

U 

U.S.A.,  Aleutian  Is.,  Shemya  I. 

52.7 

N 

174.1 

E 

5 

U.S.S.R.,  Tura,  Siberia 

64.8 

N 

101.0 

E 

6 

Finland,  Kuopio 

62.7 

N 

28.0 

TT* 

T 

Azores  Is.,  Pico  I. 

39.0 

N 

28.5 

IJ 

8 

Dutch  Guiana,  Paramaribo 

05.5 

N 

55.2. 

T.T 

9 

Equaclor,  Quito 

00.1 

S 

78.5 

TJ 

10 

Clipperton  Island 

10.3 

N 

109.2 

W 

11 

U.S.A.,  Hilo,  Hawaii 

19.8 

N 

155.0 

12 

VJake  Island 

19.7 

N 

166.2 

E 

13 

Japan,  Kagoshima 

31.7 

N 

130.6 

E 

14 

India,  Gauhati 

26.2 

N 

91.7 

IF 

U 

15 

Iran,  Sabzevar 

36.5 

N 

57.5 

E 

l6 

Libya,  Sirte 

31.7 

N 

16.4 

E 

17 

Liberia,  Roberts  Field 

06 .8 

N 

10.2 

W 

18 

Trindade  Island 

20.5 

S 

29.4 

V/ 

19 

Argentina,  Villa  Dolores 

32.0 

S 

65.1 

W 

20 

Sala-Y-Gomez  Island 

26.6 

S 

105.2 

V/ 

21 

Pukapuka  Island  ^ 

14.7 

s 

138.8 

W 

22 

Wallis  Is.,  Uvea  I. 

13.2 

s 

176.3 

W 

23 

New  Guinea,  Kikori 

07.3 

s 

144.2 

E 

2k 

Sumatra,  Palembang 

03.0 

s 

105.0 

E 

25 

Maidive  Is.,  Male 

04.2 

N 

73.3 

E 

26 

Sudan,  Juba 

04.8 

N 

31.6 

S 

27 

Southwest  Africa,  Bogenfels 

27.6 

s 

15.6 

E 

28 

So.  Sandwich  Is.,  Saunders  I. 

58.4 

s 

26.7 

W 

29 

Antarctica,  Peter  I. 

69.2 

s 

90.0 

w 

30 

So.  Pacific  Ocean,  Shoal 

41.5 

s 

148.6 

w 

31 

New  Zealand,  Queenstown 

45.0 

s 

168.2 

E 

32 

Auistralia,  Denmark 

35.0 

s 

117.3 

E 

33 

St.  Paul  Island 

38.7 

s 

77.0 

E 

3^ 

Madagascar,  Fort  Dauphin 

25.0 

s 

47.1 

E 

35 

Antarctica,  U.S.S.R.  Station 

68.0 

s 

46.4 

E 

3o 

Antarctica,  France  Station 

67.0 

s 

139.0 

E 
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II.  DISCUSSION  OF  ANALYTICAL  APPROACH 


A  brief  description  of  the  method  used  to  determine  the  propaga¬ 
tion  of  errors  over  the  net  will  show  how  the  direction  to  each 
satellite  point  is  incorporated  in  the  final  least  squares  adjustment. 
It  is  assumed  the  reader  is  familiar  with  the  theory  developed  in 
Ballistic  Research  Laboratories  Report  No.  1065.^ 

Applying  this  theory  the  following  is  assxmied. 

1.  Each  line  of  sight  belonging  to  a  specific  satellite 
point  is  imaged  at  the  center  of  the  corresponding  photograph,  i.e., 

1^  =  ly  =  0,  Xp  =  y^  =  0  and  =  0 

2.  K  =  0 

Enforcing  the  orientation  parameters  Gc  and  oo  as  are  determined 
in  practice  from  the  photographed  stars,  there  remains  only  the 
station  coordinates  X^,  Y^,  and  the  satellite  coordinates  X,  Y,  Z 

which  must  be  determined  by  a  least  squares  adjustment. 

Equation(3)^  which  expresses  the  fimctional  relation  that  exists 
for  an  individual  ray  can  now  be  expressed  by  the  corresponding 
observation  equation  of  the  form 

f(V^)  =  F(0^,X^) 

where  0^^  refers  to  the  parameters  X^,  Y^,  Z^  for  a  certain  camera 
station  i  and  Xj  refers  to  the  station  location  X,  Y,  Z  of  any  one 
individual  satellite  point  j . 

It  is  assumed  that  each  satellite  point  j  is  observed  from  either 
two  or  three  stations,  whereby  satellites  in  the  center  of  a  triangle 
are  seen  from  three  stations  while  those  on  the  center  of  station  lines 
are  observed  from  just  two  stations. 
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Since  ve  are  interested  only  in  station  positions,  for  computa¬ 
tional  convenience  all  satellite  positions  are  eliminated.  To 
accomplish  this.  It  is  necessary  at  the  time  of  elimination  that  aJ.1 
rays  to  a  specific  satellite  position  have  been  incorporated.  From 
the  geometrical  standpoint  no  one  ray  has  preference  over  another, 
therefore  each  ray  may  be  treated  independently.  It  is  necessary  to 
consider  only  one  set  of  two  equations  for  any  individual  ray  and 
apply  the  same  computation  to  all  rays  one  after  the  other. 

The  equation  which  simulates  mathematically  a  single  ray  is  the 
condition  that  object  point,  center  of  projection  and  corresponding 
image  point  are  colllnear. 

In  this  case,  therefore,  we  obtain  for  each  ray  two  equations: 

=  -  II,jiX  -  E^Y  -  F^YiZ  + 


In  matrix  notation  we  obtain  the  3 

matrix  for  the  stations 
and 

B  matrix  for  satellite  positions 

A 

From  Figure  (h) 


matrices,  namely 
X 


-^X 

-F. 

-Dy 

-F, 

F. 

X 

X 

D,, 

F. 

(X)  =  X  -  Xq 

(Y)  =  Y  -  Yq 

(z)  =  Z  -  Zq 


q  =  J  {xf  +  i'lf  +  (zf 


c  =  focal  length  =  constant 
14 


toXIl  ^ 


(X) 


(Y) 


y  {xf  +  {zf 


sin  CD  s  _ 

\/yx)^T7Y)^T7^ 


cos  a 


(Z) 

y  (x)^  +  (z)^ 


cos  O) 


y  (x)^  +  (z)^ 
y  (x)^  +  (Y)^  +  {zf 


Since  /c  =  0  the  direction  cosines  become: 


=  -  cos  Cc 


Bl  =  0 


Ag  =  -  sin  a  cos  cd  =  cos  cd 


D  =  sin  a  cos  CD 


E  =  sin  0) 


and  thus  the  computational  auxiliaries  are : 


sin  a 


cos  a  sin  cd 


cos  CC  cos  CD 
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By  substitution,  the  coefficients  of  the  observation  equations 
are  found  to  be; 


Y  (x)^  +  (z)‘ 


n  -  ^  (X)(Y) 

Y  "  2  / - = - 


y  (x)^  +  (z)^ 


Ex  =  0 


Ey  =  ^  y  (x)^  +  (z)^ 


F  =  -  E 
X  q 


(X)^  +  (Z)^ 


F  =  — 
2 


(Y)(Z) 


‘I  y  (X)*"  +  (Z)‘ 


To  form  the  normal  equations  we  must  compute  for  each  ray  the 
following  expressions : 


[dd]  =  =  \  JY)^  +  (Zf^ 


I  DE~1  =  D  E  +  D  E  = 
L  J  XX  Y  Y 


_  2/ \  \ 
c  l  AJl  I  J 


[df]  .  ^  . 


c^(X)(Zl 


[ee]  =  e|  +  e|  =  ^  (X)^  +  (z)^ 


[ef]  . 


E^F^  +  E^Fy  =  -  (Y)(Z) 

q 


[ffJ  =  F^Y  =  -^  (Y)^  +  (Y)^ 
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Since  each  ray  is  treated  independently  it  will  contribute  to 
the  normal  equation  system  the  following  submatrix: 


Satellite  Station 

AX  AY  AZ  AX^  AY^ 

3  s  s  0  0  0 


[ddJ 

[DE] 

[df] 

[dd] 

[de] 

[-df] 

[de] 

[ee| 

[ef] 

[de] 

[ee] 

[-ef] 

[df] 

[ef] 

[?f] 

[df] 

[ef] 

[ff] 

Pdd] 

[-Dd 

[df] 

[dd] 

[de] 

[df] 

[-de] 

[-EE] 

[ff] 

[de] 

[ee] 

[ff] 

[-df] 

[ef] 

[-FF] 

[df] 

[ef] 

[ff-] 

Figure  5 


Figure  5  shows  the  6x6  submatrix  for  one  ray  only.  A  satellite 
vrill  never  be  seen  from  more  them  three  stations  so  this  6x6  submatrix 
\n.ll  be  computed  a  maximum  of  three  times,  once  for  each  satellite 
ray.  If  the  satellite  is  observed  from  only  two  stations,  the  third 
6x6  submatrix  vrill  appear  as  a  null  matrix. 

For  each  individual  obser^ration  we  are  able  to  form  a  corres¬ 
ponding  partial  normal  equation  system, and  by  adding  each  additional 
ray,  the  final  normal  equation  system  can  be  accumulated  as  shovm 
in  Fig\3.res  6,  A_bCD. 

At  this  point  in  the  computational  procedure  the  satellite 
point  is  eliminated  completely. 


D  = 


B.": , 

^/o 


1  A 


-LW 


In  general 


12 


CD 

B, 

FIGURE  6A 


Figure  6A  represents  the  contribution  of  one  ray  to  the  normal 
equation  system. 
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-Bi 

-B| 

B| 

-Ba 

Be 

FIGURE  6B 


Figure  6B  represents  the  contribution  of  two  rays  to  the 
normal  equation  system. 
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12 


♦B3 

-Bg 

-B3 

-B, 

B| 

02 

-Bj 

B3 

FIGURE,  eg 


Figure  6C  represents  the  contribution  of  ail  three  rays  to  the  normal 
equation  system. 
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CD 

X 

Bx/0 

rT 

“X/D 

CD 

0 

FIGURE  6D 


The  upper  left  3x3  Bv-matrix  contains  coefficients  associated  with  the 
unknown  corrections  of  the  satellite  point,  AX,  AY,  and  AZ. 

The  lower  right  9x9  Bq" matrix  is  subdivided  into  three  block  diagonal 
matrices,  one  for  each  camera  station  that  sees  the  satellite  point,  and  is 
associated  with  the  corrections  of  the  unknown  positional  parameters  AXq, 
AYo,  and  AZ©  for  each  station. 

The  Bx/0"niatrix  contains  the  coefficients  which  represent  the  geo¬ 
metrical  fact  that  the  satellite  point  was  seen  from  stations  shown  in 
the  Bo-matrix. 
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This  elimination  cycle  reduces  the  12x12  matrix  of  Figure  6d  to 
a  9x9  D-matrix  as  shovn  by  Figure  7,  concerned  with  station  data  only. 


Sta.  i 

sta.  0 

Sta.  k 

o 

o 

o 

o 

o 

o 

o 

o 

Station  i 

^ii 

°i0 

®ik 

Station  j 

Station  k 

^ik 

^kk 

Figure  7 

In  Figure  7  the  satellite  point  has  been  numerically  eliminated, 
and  only  station  unknowns  euid  AZ^  remain  for  two  or  three 

stations  depending  on  how  many  stations  recorded  the  specific 
satellite  point. 

This  study  is  made  for  a  maximum  of  thirty-six  stations  where 
n  denotes  the  number  of  ground  stations  included  in  a  specific  net. 

The  size  of  the  reduced  normal  equation  system  will  be  a  3nx3n  matrix 
containing  unknown  ground  stations  only. 

Consequently,  each  9^9  D-matrix  is  partitioned  into  nine 
3x3  submatrices,  each  of  which  is  added  in  its  proper  position 
within  the  final  3nx3n  matrix  system  as  shoim  in  Figure  8. 

Before  the  resulting  system  is  inverted,  the  following  additional 
conditions  must  be  incorporated: 

l)  The  geometric  condition  that  sets  the  sum  of  all  adjusted 
Xq,  Yq,  Zq  coordinates  individually  equal  to  zero  in  order  to  obtain 

a  symmetrically  arranged  error  distribution  field. 
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Assume  one  satellite  point  is  observed  from  stations  3,  5,  and  7. 

The  9x9  D-matrix  is  sub-divided  and  added  into  the  3n  x  3n 
system  (n  ■  7)  as  follows: 

i  ■  3 
j  “  5 
K  *  7 


- STATION  NUMBER - 

2  3  4  5  6  7 


AXo  AYo  AZqAXo  AYqAZo  AXqAYq  AZq 


Dii 

Dij 

Dik 

“ii 

...  i 

D-.X 

1 

_ 1 

t>kk 

FIGURE  8 
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1^0=  -&0 

2)  Metric  condition:  At  least  one  scalar  must  be  introduced. 
Allowance  has  been  made  in  the  computer  program  to  incorporate  as 
many  as  four  scalars  \rith  their  respective  variances. 

If  the  distance  between  two  stations  (i  and  j)  is  enforced, 
the  corresponding  metric  conditional  equation  is  introduced: 


G'AX^  +  H-AY.  +  J-AZ.  -  G-AX.  -  H*AY„  -  J-AZ^  =  0 
Oi  Oi  Oi  Oi  Oi  Oi 


where 


G  (^0." 

10 


"  =  ‘0  Viij 

X  J 


J  .  (Zq  -  Z„  )/d 
1  ,) 


These  conditions  are  incorporated  in  the  final  normal  equation 
system  by  simply  adding  the  specific  equations  to  the  normal  equation 
system  (3nx3n)  and  at  the  same  time  restoring  the  system  to  a 
quadratic  form.  (paq,e  33 

Assuming  n=36  and  the  number  of  scalars  introduced  to  be  four, 
Figure  9  shows  the  resulting  system  after  the  geometric  and  metric 
conditions  have  been  incorporated,  A  system  of  this  size  refers  to 
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FINAL  MATRIX  TO  BE  INVERTED 

24 


108  station  coordinates  and  requires  the  Inversion  of  a  115x115 
matrix.  The  inverted  Q-matrix  refers  to  the  Cartesian  coordinate 
system  in  which  the  triangulation  was  performed  with  its  origin  at 
the  center  of  gravity  of  the  whole  system. 


Q  = 


The  diagonal  terras  can  be  used  to  compute  the  mean  errors  of 


the  station  coordinates 


with 


UV  =  \x’10^ ’I  % 


For  each  station  a  transformation  matrix  "X obtained  from  a 

coordinate  transformation  subprogram,  is  used  to  transform  the  diagonal 
terms  of  the  Q-matrix  from  the  Cartesian  coordinate  system  to  a 
geodetic  ellipsoidal  system,  in  terms  of  latitude  (0),  longitude  (X) 
and  elevation  (H). 


Appropriate  constants  are  applied  to  the  diagonal  terms  to  give 
results  in  meters. 

J 


®0  = 


m 


where  =  30.85 

where  =  30.85  cos  0  J  Q.x 
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where 


m 


H 


A  computer  program  was  written  at  3RL  in  the  Forast  programming 
language  which  enables  the  investigation  of  error  propagation  for 
various  size  nets  up  to  a  maximum  of  36  stations,  allowing  the  rramber 
of  scalars  to  vary  from  one  to  four. 


An  error  analysis  of  various  station  schemes  using  hypothetical 
data  is  presented.  The  scheme  which  is  thought  to  be  most  nearly 
ideal  consists  of  36  stations  forming  68  triangles  with  102  sides, 
4000  to  4500  km  in  length,  using  a  satellite  altitude  of  3600  km. 

In  this  scheme  it  was  assumed  that  all  satellite  positions  over  the 
center  of  each  of  the  66  triangles  were  observed  simultaneously  from 
three  ground  stations  and  satellite  positions  over  lines  were 
observed  from  two  stations.  A  total  of  I70  satellite  positions  and 
4o8  rays  were  incorporated  in  the  least  squares  adjustment  of  the 
36  station  net.  Four  specific  nets  are  considered: 


No .  of 
stations 

in  net 

liO.  of 
triangles 
formed 

— Wo.  of — 
lines 

No .  of 
satellite 
positions 
observed 

NO.  or  rays 
incorporated 
in  least  square 
solution 

36 

68 

102 

170 

106 

35 

62 

96 

l64 

390 

3U 

56 

90 

156 

372 

26 

36 

66 

122 

282 

T  c?  4- ri  4- 

•JUC3.0XWAX  11  u  oiac.  ij  o  x  v-^iio 


CXX  1  CLll^^ 


the  entire  earth. 


2.  In  the  35  station  net  Russia  is  omitted. 


3.  In  the  34  station  net  Russia  and  the  South  Pacific  are 
omitted . 


4.  In  the  26  station  net  the  stations  in  the  polar  regions 
are  omitted  leaving  a  28  station  equatorial  net  with 
two  rows  of  triangles  straddling  the  equator. 
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The  propagation  of  the  errors  in  the  station  coordinates  of  the 
aforementioned  triangulation  nets  are  presented  in  graphical  form. 
These  results  have  been  obtained  by  inverting  the  final  normal 
equation  system.  The  diagonal  terms  of  the  Inverted  matrix  give 
quantitative  information  about  the  propagation  of  errors  in  the 
unknowns  of  the  solution  and  the  non-diagonal  terms  furnish  the 
corresponding  covariances.^ 

In  all  schemes  the  effect  of  a  specific  number  (one  to  four)  of 
scalars  and  various  accuracies  of  these  ground  established  baselines 
on  the  propagation  of  errors  over  the  net  were  investigated . 

The  stations  between  which  distance  measurements  are  assumed 
are  listed  below  and  are  indicated  by  scalar  number  on  the  vertical 
scale  of  Figures  10-49. 


Scalar  No. 

Baseline 

Stations  between  which 
distance  measurements 
are  assumed 

1 

North  America 

2  -  3 

2 

Australia 

23  -  32 

3 

Europe 

6  -  l6 

3* 

Iran 

l4  -  15 

4 

South  America 

8  -  19 

*  For  the  twenty-eight  station  net  the  third  scalar  between 
stations  6  and  l6  was  replaced  by  a  scalar  in  Iran  between 
stations  l4  and  15* 
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III.  DISCUSSION  OF  RESULTS 


In  order  to  introdiice  the  various  scalars  with  different 

variances  it  is  necessarj^  to  compute  the  corresponding  weighting 
7 

factors. 

Assxjming  a  standard  error  of  unit  weight  on  the  plate  measure¬ 
ments  of  t  p  microns  we  have 


K 


P  =  1  =  - 

^  2  .  -12 
p  "10 


2  -12 

or  K  =  p  *10 


Correspondingly  the  weight  for  any  one  introduced  scalar,  if 
its  specific  standard  error  is  denoted  by  t  m  (meters),  is 


2  -12 
K  ^  [I  -10 


m  m 


K  = 


Stipulating  a  realistic  value  for  p  of  t  0.75  microns, 

-12 

0.5 *10  auid  consequently  we  obtain  the  follovring  values  for  m. 


A 

m 

U  -It 

10  -10  ^ 

assumed  to 

be  zero 

-Ik 

1-10 

±  T  m 

-it 

0.1-10 

±  20  m 

For  convenience  the  problem  is  scaled  by  10  resulting  in  the 

I). 

P  values  of  10  ,  1,  and  0,1,  respectively.  All  results  are  related 

to  a  focal  length  of  the  photograirmetric  camera,  c  =  0,3l6m. 

Figures  10,  11,  12  are  concerned  \T±th  the  36  station  net.  The 
accuracies  are  shoim  which  can  be  obtained  for  the  determination 
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of  station  coordinates.  The  assumed  standard  deviations  on  the 
baseline  measurements,  as  outlined  above,  are  assumed  to  be  zero, 
t  7  meters  and  ±  20  meters,  respectively.  The  beneficial  influence 
of  one,  two,  three  and  four  scalars  on  the  accuracy  of  the  determina¬ 
tion  of  the  station  coordinates  is  demonstrated.  The  geometry  of 
the  triangulation  is  strengthened  with  each  additioneJ.  scalar  and, 
consequently,  the  resiilts  show  an  improvement  in  the  mean  error  of 
the  station  coordinates. 

Additional  scalars,  even  when  measured  with  extreme  accuracies, 
make  little  improvement  on  latitude  and  longitude  determination.  In 
fact,  for  all  practical  purposes  one  scalar  is  nearly  as  good  as 
four.  However,  the  height  determination  is  decidedly  improved  with 
each  additional  scalar  and  with  an  increase  in  accuracy  of  scalar 
measurements.  The  ratio  of  the  improvement  in  the  mean  error  of 
station  coordinates  becomes  smaller  with  the  introduction  of  each 
additional  scalar.  The  greatest  relative  improvement  exists  when 
two  scalars  are  used  to  constrain  the  triangulation,  and  the  addition 
of  a  third  scalar  will  yield  significant  improvement.  However, 
whether  or  not  much  is  gained  from  four  scalars  over  three  is 
questionable.  This  fact  is  significant  when  considering  the  cost  of 
obtaining  baseline  measurements  over  long  distances.  The  influence 
of  the  measuring  errors  as  they  are  propagated  over  the  net  can  best 
be  seen  by  examining  Figure  13  which  is  a  combined  plot  of 
Figures  10,  11,  12.  Measuring  errors  of  1  7  meters  using  one  baseline 
give  about  the  same  accuracy  as  four  baselines  with  measuring  errors 
of  t  20  meters.  Similiarly,  for  one  scalar  with  an  assumed  standard 
deviation  of  t  20  meters  the  errors  propagated  over  the  net  are 
twice  as  large  as  those  obtained  with  the  corresponding  baseline 
measurements  of  t  7  meters  accuracy.  The  table  that  follows 
illustrates  the  situation.  The  fact  should  not  be  overlooked  that 
only  an  insignificant  improvement  is  obtained  between  assumed  zero 
measurements  and  those  with  an  assumed  accuracy  of  t  7  nieters . 
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-  values  for 

■if 

H-determlnation 

No.  of 

Assumed  standard  deviation 

of  scalar 

scalars 

1 20 

-  7 

0 

1 

5.8 

3.1 

2.7 

2 

k.2 

2.k 

2.1 

3 

3.6 

2.1 

1.9 

k 

3.1 

1.9 

1.8 

* 

Portions 

of 

these  results  were  presented 

.  8 
in  . 

The  possibility  of  being  unable  to  establish  certain  station 
locations  led  to  the  investigation  of  several  other  nets.  The  same 
identical  runs  that  were  computed  for  the  36  station  net  were  computed 
for  the  following  nets: 

a)  35  station  net  without  Russia  (station  No.  5  omitted) 

Figures  14-17. 

b)  3^  station  net  without  Russia  and  the  South  Pacific  (stations 
No.  5  and  30  omitted)  Figures  18-21. 

c)  26  station  net  without  stations  near  the  polar  regions. 
Figures  22-25. 

For  each  station  not  incorporated  in  the  least  squares  adjustment 
the  triangulation  is  vreakened;  this  accounts  for  the  more  erratic 
appearance  of  each  set  of  Figures  for  the  35>  34  and  28  station  nets, 
respectively. 

In  the  least  squares  adjustment  for  the  35  station  net  six 
triangles  in  the  vicinity  of  Russia  are  omitted.  The  irregularity 
of  Figures  14-17  for  all  three  coordinates,  latitude,  longitude  and 
height  in  this  particular  area  of  the  net  is  appareni  ,  particularly 
at  stations  which  are  close  to  the  latitude  aid  longitude  of  the 
omitted  station  No.  5*  Station  6,  which  is  about  the  same  latitude 
as  station  5^  has  the  largest  mean  error  in  latitude  over  the  net 
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while  stations  13  and  lU,  v/hich  have  longitudes  near  that  of  station  5> 
have  the  largest  mean  er:>'’or  in  longitude  over  the  net.  Furthermore, 
the  height  determinations  have  the  largest  mean  errors  occurring  on 
stations  6  and  l4  which  are  closest  in  both  latitude  and  longitude 
to  station  p.  On  the  other  hand  the  loss  in  accuracy  is  relatively 
small,  on  accoimt  of  the  basic  favorable  error  propagation,  which 
is  unique  for  this  tyi^e  of  three-dimensional  satellite  triangulation. 

In  Figures  26-49  the  same  information  is  presented  in  a 
different  way,  by  comparing  the  accuracies  obtainable  for  the  four 
different  net  sizes  (3t>,  35  34  and  26  stations)  in  terms  of 
latitude,  longitude  and  height. 
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